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Abstract 

We consider systems, which conserve the particle number and are described 
by Schrodinger equations containing complex nonlinearities. In the case of 
canonical systems, we study their main symmetries and conservation laws. 
We introduce a Cole-Hopf like transformation both for canonical and non- 
canonical systems, which changes the evolution equation into another one 
containing purely real nonlinearities, and reduces the continuity equation to 
the standard form of the linear theory. This approach allows us to treat, in 
a unifying scheme, a wide variety of canonical and noncanonical nonlinear 
systems, some of them already known in the literature. 
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I. INTRODUCTION 



Over the last few decades many nonlinear Schrodinger equations (NLSEs) have been 
proposed in order to test the fundamental postulates of quantum mechanics, like for in- 
stance, the Bialynicki-Birula and Mycielski equation ffl], the Kostin equation [0], the Gisin 
equation |1 among many others M. In Ref. [H a wide class of NLSEs for finite-dimensional 
quantum systems was selected in order to preserve the homogeneity principle of the origi- 
nal Schrodinger equation, with the superposition principle being destroyed by the nonlinear 
terms. 

Many of the NLSEs proposed in literature contain complex nonlinearities. For instance, 
the Doebner-Goldin (DG) equations |6|-[8[ were introduced as the most general class of 
Schrodinger equations, compatible with the Fokker-Planck equation for the probability den- 
sity p = \if)\ 2 namely d p/dt + V • j = D A p, being j the standard quantum current and 
D a positive diffusion coefficient. The importance of this class of evolution equations is that 
it is founded on the grounds of the group theory: the nonlinear terms were derived from the 
representation analysis of the Diff(R 3 ) group which was proposed as a universal quantum 
kinematical group |§ . 

In addition, a large number of NLSEs with complex nonlinearities have been proposed 
in order to describe some phenomenologies in condensed matter physics. For instance, in 
Ref. |10| a NLSE with a nonlinearity of the type a\ \ip\ 2 ip + (12 |V , | 4 V ; + i a 3 d x (\ip\ 2 ip) + 
(04 + i as) d x \ijj\ 2 ip is introduced to describe a single mode wave propagation in a Kerr 
dielectric guide. Another example is given by the generalized Ginsburg-Landau equation 



appeared in Ref. |TTJ. This equation contains the nonlinearity: a\ \ip\ ip+i ct2 ip + i d xx ip + 
i 04 |^| 2 ip 0] which takes into account pumping and dumping effects of the nonlinear media 
and can be used to describe dynamical modes of plasma physics, hydrodynamics and also, 
solitons in optical fibers (Ref. |13[ and references therein). Finally, complex nonlinearities 
in Schrodinger equations are used to describe propagation of high power optical pulses in 



ultrashort soliton communication systems [0,|T^] , incoherent solitons, ||T6|,|TT1 , multi-channel 
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bit-parallel-wavelength optical fiber networks ||18|| , among others. 

In this paper we consider the most general class of NLSEs conserving the quantity iV = 

f\*fj\ 2 d n x: 

ih^ = -^-^ + {w + iW)ii , (l.i) 

at 2m 

where the real W and imaginary W parts of the complex nonlinearity are smooth functions 
of the fields ip, ip* and their spatial derivatives of any order. When ip is written in polar 
representation ip = p 1//2 exp(i S/h), Eq. (|1 . 1| ) is split into two nonlinear partial differential 
equations for the real fields p and S: 

»£ + V.(™, + fUo, (1.2, 



dt \ m 

OS (VS) 2 



dt 2m 

where: 

h 2 A 



+ W + U g = 0, (1.3) 



"' = -4^Vzf' (L4) 



is the quantum potential |19| and the real functional F is related to W through 



W = Av-F, (1.5) 
2p 



as the particle number conservation requires. It is easy to recognize that Eq. ( |1.2| ) is a non- 
linear continuity equation, which involves only the term W, while Eq. (O) is a nonlinear 
Hamilton-,] acobi like equation involving only the term W . 



In the Calogero picture ||20|| , the system ( |1.2| ), (O) is C-integrable if there exists a trans 



formation of the dependent or/and independent variables: t —>■ T, x —>■ X, p —>■ R, S —>■ S 
which transforms Eqs. ( p..2|) , ( |1.3|) into: 
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V and U q being the gradient and the quantum potential in the new variables. Eqs. 
( P-.6Q , ( |1.7|) constitute the well known hydrodynamic representation of the standard linear 
Schrodinger equation. 

The principal aim of this paper is to introduce a nonlinear transformation for the field 
S: S —>■ S, in order to reduce Eq. ( |1 . 2| ) to the standard form of the linear theory ( |1.6|) . As 
a consequence of this transformation, the evolution equation ( |1 . 1| ) transforms into another 
one containing a purely real nonlinearity. Moreover, the current, when expressed in terms 
of the new field <fi = p 1 ^ 2 exp(iS/h), reduces to the standard bilinear form of the linear 
Schrodinger theory. 

The paper is organized as follows. In Sect. II, we introduce a general class of (n + 1) 
canonical NLSEs, invariant over the action of the U(l) group. In Sect. Ill, starting from 
the Noether theorem, we consider the main symmetries and related conserved quantities 
of the canonical system. In Sect. IV, we introduce a Cole-Hopf like transformation which 
eliminates the imaginary part of the nonlinearity in the evolution equation, while in Sect. 

V, the same transformation is considered in the case of a noncanonical systems. In Sect. 

VI, in the framework of the approach developed in the previous sections, we treat, in a 
unifying context, some NLSEs already known in the literature, in order to show that all the 
transformations introduced by the various authors to study these equations can be obtained 
as particular cases of the transformation here proposed. Finally, some conclusions and 
remarks are reported in Sect. VII. 



II. THE CANONICAL MODEL 

Let us consider the class of canonical NLSEs described by the Lagrangian density: 

where V = (d±, • • • , d n ) is the n-dimensional gradient operator. The last term in the r.h.s. 
of Eq. ( |2.1|) is the nonlinear potential which we assume to be a real smooth function of 
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the fields ip and ip* an d their spatial derivatives. Here and in the following, we use the 
notation U[a] to indicate the dependence of U on the field a and its spatial derivative of any 



order. We deal with dynamical systems described by Eq. ( |2.1| ) which are invariant under 
transformations belonging to U(l) group. As we will show in the next section, this condition 
imposes a constraint on the form of the nonlinear potential U. 
We start from the action 



A = Cd n xdt 



(2.2) 



and observe that the evolution equation of the field xp is given by: 

8 A 



S ip* 



. 



(2-3) 



The functional derivative is defined through |pTj| : 



5_ 

5 a 



Jg[a}d n x= E(-i) fc ^ 



[fe=0] 



dG[a] 
d(p,a) 



(2.4) 



with T> t = d k /(dx l i 1 ■ ■ -x 1 ™) and J2[k=o] = 12kLo SjV The sum J2i k is over the multi-index 
I k = i 2 , ■ • ■ , i n ) where 1 < p < n, < i p < k and = k. 
Eq. (|2.3|) assumes the form: 



(2.5) 



which, after performing the functional derivatives, transforms to the following NLSE: 

dip h 2 . , 5 



i h 



AiP + - — / U[ip*, ip]d n xdt 
ow* J 



(2.6) 



where A = d 1 
form: 



dt 2m T ' Sip* 
d\ is the Laplacian operator. Eq. (|2.6|) can finally be written in the 



(2.7) 



where the complex nonlinearity 



A[p, S} = W[p, S}+iW[ P , S] , 



(2.8) 



has real W[p, S] and imaginary W[p, S] part defined through: 



W\p,S\ = -=- [ U[p, S]d n xdt , 
o p J 

W[p, S] = f U[pj S ]d n xdt 

2 po b J 



(2.9) 
(2.10) 



p and S being the hydrodynamic fields related with the wavefunction ip through [19,^2 



ip(x, t) = p 1/2 (x, t) exp 



tS(x, t) 
n 



[2.11) 



III. SYMMETRIES 



In this section we study the main symmetries and conserved quantities of the system 
described by the Lagrangian (|2.1|) . 

Let us consider the U(l) invariance condition. The variation 5 e ip = ieip, with e an 
infinitesimal real parameter, implies the following variation on the action: 





5 e A = -eh f jr^U[p, S] d n x dt . 

J O 



(3.1) 



Taking into account the Noether theorem [£3(] we can also write the variation 5 e A in the 
form 



5 e A= -eh J d^j^ip*, ip) d n xdt 



(3.2) 



By comparing Eq. ( [3.1|) and (|3.2j ) we obtain: 



dp „ . dU 

dt J dS 



(3.3) 



with p = jo- Eq. ( |3.3| ) is not a continuity equation because the Lagrangian ( |2.1| ), for a 
general nonlinear potential U[p, S], is not [/(l)-invariant. In appendix B we show that 
£/(l)-symmetry can be restored if one assumes that the nonlinear potential U[p, S] depends 
on S only through its spatial derivative, modulo a total derivative term, which does not 
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change the dynamics of the system (null Lagrangian). As a consequence, the r.h.s of Eq. 
( p.3| ) vanishes and it becomes a continuity equation for the conserved density p. Thus, 
the [/(l)-invariance limits the class of nonlinear potentials appearing in Eq. ( |2.1| ). In the 
following, we consider only £7(l)-invariant systems, where the functional U[p, S] depends on 
S trough its spatial derivative. For this class of systems, Eq. (|3.3|) becomes: 



!? + V.J = 0, (3.4) 
The conserved quantity associated to the continuity equation (|3.4j) is: 

N = [ pd n x . (3.5) 



The expression of j is obtained in appendix A, and is given by: 

au\p,s}' 



m [k=0] /, 

where 



d(V iJt S) 



(3.6) 



n 



fi k+1 = n ~Yl <W , (3.7) 

for l h = (rrii, ni2, • • • , nT>n) with 1 < r < n, < m r < n and J2r m r = k + 1. 

Note that the expression ( p.6|) of j, can also be written (see appendix A) in the form: 



ji = —p+—^— I U[p, S]d n xdt . (3.8) 



m r 6(diS) 

Eq. ( p.8|) can be obtained starting directly from the equation (p.7|) after adopting the 
hypothesis that U[p, S] depends on the field S only through its spatial derivatives as required 
from the U(l) symmetry. 

In the following we consider the main space-time symmetries of the Lagrangian fl2.1| ). 



We note that U[p, S] depends on the variables x and t only trough the fields p and S, thus 
the system is invariant over space-time translations. From Noether theorem we have: 

^ + V-T, = 0, (3.9) 



where 7^ = T ^; (T ^)i = T ifl with fi = 0, • • • , 3. The components of the energy-momentum 
tensor T^ v (see appendix A) are given by: 

n 2 



00 



2m 

h 



\V^\ 2 + U[p, S] 



(3.10) 
(3.11) 



T, 



,0 



2 m 



{diip* d t ip - diip d t ijj* 



T 



'J 



E E(-iK^ 9 <K d1 7 [p ' s \ 

[fc=0][p=0] " ( " [ d \PrJ k P) 

Ti 2 

(diil)* djifj + djip* diip) + 5ijC 



dU[p, S] 

9 & Ik s) 



d t S 



(3.12) 



2 m 



E E(-i) p <2\R 

[k=0] [p=0] 



dU\p, S] 
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From Eq. ( |3.1(JD and ( |3.11| ) we obtain the conserved quantities: 

ti 2 



9 K*) 



(3.13) 



E 
P 



2 m 
h 



\ViP\ 2 + U[p, S] 



d n x , 



(3.14) 
(3.15) 



which are respectively the total energy and the linear momentum of the system. We can see 
that U[p, S] modifies the expression of the energy while the momentum maintains the form 
of the linear theory. From Eqs. ( |3.12| ) and (|3.13|) we see that the presence of U[p, S] also 
modifies the expressions of the fluxes associated to E and P. 

Note that if the energy-momentum tensor is symmetric in the spatial indices = Tj i: 
the potential U[p, S] is invariant over the action of the orthogonal group SO(n). In this 
way we can define n (n — l)/2 conserved quantities: 



a^-,a n - 2 _ e ax,-,an_ a ,ij j x . Tq .^ X 



(3.16) 



where e ai '"'^ is the n-rank totally antisymmetric tensor defined as e 1 ''"' 1 = 1. For n = 3 we 
recognize the well known conserved components of the angular momentum. 

Finally, we look at the Galilei invariance. We recall that if the system admits this 
symmetry, the corresponding generator: 
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G = Pt-mNxr 



(3.17) 



is conserved, namely: 



d_G 

~dt 



(3.18) 



where the linear momentum is given by: 



P = J p VSd n x 



(3.19) 



while the mass center vector is defined as 



wJ px<rx 



(3.20) 



We consider now the Ehrenfest relation: 



dx c 

at 



(3.21) 



which can be obtained from Eq.s (|3.4j ) and (|3.20|) . We note the formal similarity with the 
corresponding relation of the linear theory. Here the expression of j is given by Eq. ( |3.6| ) and 
depends on the form of the nonlinear potential U[p, S]. Taking into account the conservation 
of P, and after assuming uniform conditions, from Eq. (|3.17 ) we obtain: 



dG l 
~~dt 



-m e (-if^r 1 /^ 



[fc=0] 



aujp, s] 



(3.22) 



Eq. ( |3.22| ) shows that the presence of U[p, S] breaks the Galilei invariance, which can be 
restored if W = as we can verify easily by using Eq. ( |2.10p . 



IV. A COLE-HOPF LIKE TRANSFORMATION 



Let us introduce a unitary transformation of the field ip'- 



i){x, t) -> <p(x, t) = U[p, S] if>(x, t) 



(4.1) 



with U* =U 1 so that: 



M 2 = \4>\ 2 = p 



(4.2) 



The functional U is chosen to eliminate the imaginary part of the NLSE ( |2.7|) and, at the 
same time, to transform the current j, given by Eq. ( |3.6| ), into another current j —>■ J 
having the canonical form: 



V5 

J = p , 

m 



with S being the phase of the new field 4>: 

4>(x, t) = p 1 ^ 2 (x, t) exp 
We write U[p, S] as follows: 



(4.3) 



-S{x, t) 



(4.4) 



U[p, S] = exp [-^<r\p, S] ) , 



(4.5) 



and observe that the generator a[p, S] is a real functional which allows us to calculate S 
starting from S: 



S = S + a[p, S] . 

The generator a[p, S] can be obtained easily by combining Eq.s ( |3.8|) , ( f4.3|) and 

m 5 



(4.6) 



di<r[p, S] 



U[p, S]d n xdt . 



(4.7) 



Eq. (|4.7|) imposes a condition on the form of the nonlinear potential which can be obtained 
using the relation: dy a = dji a: 



5 



0; 



s 



J U[p, S}d n xdt = . 



(4.8) 



Condition ( |4.8j ) selects the potentials U[p, S] and the nonlinear systems where we can per- 
form the transformation (|4.1|) . In the case of one-dimensional systems the transformation 
Q4.1D is always accomplished. 

It is now easy to verify that the transformation ( O ) reduces the evolution equation ( ^.7| ) 
to the following NLSE: 
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which now contains only the real nonlinearity W[p, S] given by : 



W[p, S] = W + 




J -Vex da 



(4.10) 



2m 



dt ' 



where W = W[p, S[p, <S]]. The phase S appears in Eq. ( |Q|) only trough its spatial 
derivatives, consequently the arbitrary integration constant, deriving from the definition of 
U, does not produce effects and can be posed equal to zero. Note that W depends implicitly 



where the current J now takes the standard expression of the linear quantum mechanics 
given by Eq. (|4.3|) . 

In conclusion, we have introduced a nonlinear and nonlocal transformation which makes 



equation (|3.3| ) to the bilinear standard form. The price that we pay is that Eq. (|4.9| ) is not 
generally canonical because the transformation (|4.1|) is itself not canonical. 

We briefly discuss the conditions under which the system, described by Eq. (f4.9|) , be- 
comes a canonical one. The canonicity of the system implies the existence of a nonlinear 
potential U from which we can derive the nonlinearity of Eq. ( |4.9|) . 

We observe that the absence of the imaginary part W in the nonlinearity of Eq. (|4.9|) re- 
quires that U depends only on the field p and its spatial derivatives. Consequently, W is a 
functional of the field p linked with U (p) through the relation 



on the field S. In fact, Eq. ( |4.6|) defines S as a function of p and S. 
From Eq. ( f4.9|) we can obtain the following continuity equation: 




(4.11) 



real the complex nonlinearity in Eq. ( |2.7| ) and at the same time reduces the continuity 




(4.12) 



which after performing the functional derivative assumes the form: 



w\p] = E (-1)*? 



[fe=0] 



" dU[p] 
k [d(V Ik p) 



(4.13) 
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In section VI we will consider few particular systems where the condition (|4.13|) is satisfied 
and their canonicity is preserved. 



V. NONCANONICAL SYSTEMS 

In this section we consider the transformation introduced previously and study its ap- 
plicability in the case of noncanonical systems. 

For noncanonical systems, the evolution equation is given by Eq. (LI]), where W[p, S] 
is now an arbitrary functional, while W[p, S] assumes the form: 

W[p,S] = ^-V-F[p,S], (5.1) 
z p 

enforced by the conservation of N = J pd n x, with F[p, S] an arbitrary functional, and the 
current is given by: 

j = —p-F\p,S\. (5.2) 
m 



It is easy to verify that the transformation ( f4.1| ) eliminates the imaginary part of the 
nonlinearity in the motion equation, which transforms again into Eqs. (t4.9|) - (|4.10|) . The 
generator a of the transformation ip — > is related to F through: 

Tfl 

V a[p, S] = - F[p, S] , (5.3) 
P 

while the condition 

F 

V x — = , (5.4) 
P 

permits the definition of the transformation in any n > 1 spatial dimension. This condition 



constrains only the form of W, differently to the canonical case, where the condition Q4.8| ) 
constrains the form of the nonlinear potential U and consequently, both W and W in the 
motion equation. 
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VI. EXAMPLES 



In this section we consider some equations already known in literature in the framework 
of the approach developed here. We show how the nonlinear transformations proposed 
to study the various NLSEs, can be obtained in a unified way as particular cases of the 
transformation given by Eq. 

i) As first trivial example we consider the canonical NLSE introduced in Ref. |[L6|| : 

ih^ = -—A^j+ (p p -^.vs)ij + - (a- Vlogp) V , (6.1) 



dt 2m V h J 2 

with 3 and a real and constant arbitrary parameters. This equation can be derived from a 
Lagrangian containing the following nonlinear potential: 

U[p,S} = ^p 2 - P -a-VS . (6.2) 

The transformation with generator a, given by: 

m 

a = — — ax, (6.3) 

produces the new canonical evolution equation: 

^ deb K 2 . , „ , ma 2 , . lS 



with associated nonlinear potential: 



8 2 ma. 2 
2 P + ^2h 

Eq. (|6.4j) can be reduced to the cubic NLSE by changing the phase S ^ S — 3ma 2 t/2h 



ii) Let us consider as second example of canonical NLSE the Chen-Lee-Liu equation \ 2% : 



dt 2m d x 2 h dx P 2 d x 
with a a real coupling constant. The nonlinear potential associated to this equation is: 
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tt\ rn a dS 2 



while the transformation with generator: 



(6.7) 



am . 
a = — -n- I pax 



2h 



reduces Eq. ( |6.6| ) to the following noncanonical NLSE: 



d(j) d 2 4> a ( d S ^ 3am \ 
dt dx 2 h \dx 8h P ) P 



(6.8) 



(6.9) 



We note that the transformation with generator a given by Eq. ( |6.8| ) is a particular case of 



the Kundu transformation introduced in Ref. [53 



iii) As third example we consider the canonical NLSE introduced in Ref.s |26| . |27 

h 2 d 2 ^ x fx 



dip 
th Jt 



dS\ , hX dp , 

o a 2 \qP + 7T } P^ + t Tr~ ■ 

2m ox 1 m V 8 ox I 2m ox 



The associated potential is given by: 



U[p, S] 



3 A 2 



i m 



p - 



X dS 



2m dx 



while the generator a assumes the form: 



a 



\\pdx 



The evolution equation for the field <fi becomes: 



i h 



90 
~dt 



h 2 <9 2 , 



A dS 



2m dx 2 m dx 



(6.10) 



(6.11) 



(6.12) 



(6.13) 



iv) As fourth example we consider the canonical NLSE recently introduced in Ref. [ f28| , f29H 



dip 



h 2 d 2 ip 



2m dx 2 m P \ d x , 



ip — i 



kTi d 



dS s 



2m p dx 

The nonlinear potential associated to Eq. (|6.14| ) is given by: 



dx 



(6.14) 



U[p, S] 



k ( dS 
2m \ dx 



(6.15) 
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Although Eq. (|6.14|) can be generalized in any spatial dimension, it is easy to verify that 
the condition ( |4.8| ) is not satisfied in this case, thus we can apply the transformation only 
to the one-dimensional case. We perform the transformation generated by: 



a = k p — — ax 

ox 



and Eq. (|6.14|) transforms into: 

dt 2 m dx 2 ml + Kp\dx 



h 2 d 2 logp 
K 1 — P ' 



(6.16) 



(6.17) 



4 m dx 2 

v) As last example of canonical system we consider the sub-class of DG equations given by 



at 2m m 



with associated a nonlinear potential: 



Ap 1 /Vp" 

T ~ 2 \ T J 



U[p, S] = apAS + P 



h 2 (Vp) 5 

m p 



The generator a now is: 



ip + ta- — %p 
2 p 



(6.18) 



(6.19) 



a = —ma logp 



(6.20) 



while the evolution equation for the field becomes: 



% n — — 



dt 2m 



A0 + 7 



Ap _ 1 / Vp " 



(6.21) 



with 7 = ma 2 — 2 (3fi 2 /m. This equation is again canonical with nonlinear potential: 



U[ P ] 



7 (Vp) ; 



2 P 



(6.22) 



and can be linearized performing the rescaling S — > S yj 2 m 7 /Ti 2 — 1 , as noted in Ref. [30 
vi) The most general class of DG equations is non canonical and takes the form: 



r\ I j- 2 5 fz 

th-^ = -—A^ + hD> >Z c i R i[P, S}^ + z-DR 2 [p, S]^ 

~~ '' ' 1 = 1 ~~ 



(6.23) 
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where R x = V-j/p, R 2 = Ap/p, R 3 = (j/p) 2 , R 4 = j-Vp/p 2 , R 5 = (Vp/p) 2 . Note that Eq. 
( |6.18| ) is obtained when D — a, Cx — — c 4 = ma/hD', c 3 = and c 2 = — 2 c 5 = — 2 (3 Ti/mD' . 
The same generator a given by Eq. ( |6.20| ) defines the transformation ip — > <p reducing the 
evolution equation to: 



dcf) h 2 5 

= -7^A<f) + Y,CiRi[p, S](j) , 

i=l 



(6.24) 



dt 2m 

where now the coefficients are given by C\ = TiD' c x — mD, c 2 = hD' c 2 , £3 = HD' C3, £4 = 
h D' C4 + m D, £5 = hD' + mD 2 /2. Note that in Ref. a nonlinear transformation was 
introduced with generator: 



a DG = ^-\ogp +} \[X(t)-l}S + e(t, x) 



(6.25) 



which produces a group of transformations mapping the DG equation into itself. We observe 
that, after posing in Eq. ( |6.25| ) 9(t,x) = 0, A(t) = 1 and j(t) = 2mf3/%, we obtain the 
generator a given by Eq. ( |6.20| ). 

vii) As a second example of noncanonical system we consider the equation: 

dip 



i h 



dt 



h 2 d 2 ifj . ,,*dij) , dip\ , 
2m ar ox ox 



(6.26) 



with q a real parameter. We note that for q = 1/2, Eq. fl6.26j ) reduces to the Kaup-Newell 
equation [31], while for q = we obtain the Chen-Lee-Liu equation ( |6.6|) . Finally for q = — 1 
the nonlinearity in Eq. (|6.26| ) becomes purely real and the equation coincides with Eq. (|6.13|) 
obtained previously. The generator: 



a 



m a 



(q + 1) J pdx 



(6.27) 



defines a transformation which reduces Eq. (|6.26| ) to: 



i h 



h 2 d 2 <\> a 



dS la 2 
^- q) dx- + 8 m h i3 - 2q - 5q)p 



P0 . 



dt 2m dx 2 h 
In the particular case q = 1 Eq. ( |6.28|) becomes canonical with nonlinear potential: 



(6.28) 



U[ P ] 



m a 
~6~F 
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P 



(6.29) 



viii) As last example we consider the Eckaus equation |32| which is a noncanonical NLSE: 



dip h 2 d 2 ip dp 2 

th ^t = -2^d^ + ia d-x^ + P p (6 - 30) 

The generator: 

ma f 

a = — J pax , (6.31) 

defines the transformation reducing Eq. (|6.30|) to the well-known quintic NLSE: 



dt 2m dx 2 \2h 
Eq. ( |6.32|) is a canonical one with nonlinear potential: 



U[P\ = ~ + P) P 3 , (6-33) 



3 V 2h 

In the particular case where (3 = —ma 2 /2h 2 the transformation with generator given by 
Eq. flOID linearizes Eq. (|O0| ). 



VII. CONCLUSION 

In this paper we have consider a class of canonical NLSEs containing complex nonlin- 
earities and describing [/(l)-invariant systems. For these systems we study the symmetries 
and the conserved quantities associated to roto-translations and Galilei invariance. 
Subsequently, we introduce a Cole-Hopf like transformation tp — > Utp, which preserves the 
f/(l)-invariance of the system and reduces the complex nonlinearity into a real one so that 
the continuity equation assumes the standard bilinear form. This transformation generally 
does not conserve the canonicity of the system. Extension to noncanonical equations is also 
studied. 

The general Cole-Hopf-like transformation introduced here, allows us to deal in a unifying 
scheme several NLSEs already known in literature, obtaining, in this way, the transforma- 
tions introduced by various authors. 
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APPENDIX A. 



In this Appendix we recover, by using the Noether theorem 
associated to a given symmetry. 
Let us consider the action: 



the continuity equation 



A 



C d n x dt 



(A.l) 



with Lagrangian: 



where 



C 



L = £ 4- C 



2 V dt * dt j 2m 



(A.2) 



(A.3) 



the standard Lagrangian density of the linear Schrodinger theory while 



-U\p, S] 



(A.4) 



is a real scalar functional depending on the hydrodynamic fields p, S and their spatial 
derivatives. The evolution equation for the field ip is given by: 

OA 



dip* 



. 



(A.5) 



Taking the functional derivatives Eq. ( |A.5| ) becomes: 



d dC T 



d dC T 



dtp* dt d (d t ip*) dxi d (diifj* 



E (-i)X 

[fc=0] 



dC, 



4 



+z — y (-i) k v r 

z ^ [k=o] 



dC, 



x/j = . 



(A.6) 



We compute the variation 5 e A generated by a one-parameter transformation group. For 
simplicity, we assume that the symmetry group acts only on the internal degrees of freedom 
of the system. The contributions to 5 e A, given by the variation of the volume element 
d n xdt, when the symmetry involves the space-time variables, are well known and can be 
added successively. Thus, we have: 



5 e A 



Sip 

dC, 

dip 



5 e ip + 



+ 



dC, 



Sip* 

d{d t tl>) 
dC 



Sp 



ss 



d n xdt 



dC 



dip* d (d t ip*) 

<9A 



d (diip*) 



+ E 

[*=0] 



d A 



5 e (V Ik S) 



d n xdt 



(A.7) 



with Z> 7 = d k / (dx 1 * ■ ■ ■ x 1 ™) and the Einstein convention for the repeated indices is assumed. 
In Eq. ( |A.7| ) we have posed J2\k=o] = ^T=o where the second sum is performed on the 
multi-index I*. = i 2 , • • • , i n ) with < i p < k, ip — k. If we use the identity: 



dC 



d_ 

da 



d{d a <P) 

with a = t, i, Eq. ( |A.7| ) becomes: 

dC, 



dC 



d{d a <j>) 



5e<P 



d_ 

da 



dC 



d(d a 0) 



(A.8) 



6 e A 



dC Lg ^ + d_ 
dtp dt 



+ 



d 

dxi 
d 

dx,. 



dC T 



d{d t iP) 
dC 



8 e iP 



d_ 
dt 

d 



+ E 

[fc=0] 



d(diip) 

dC L s 
d (diip*) 
dC, 



dtp* ^ dt 



dC 



d 

dx,. 



dC T 



d (diip*) 



d (d t ip*) 

6 e r 



d 


\ } 


~ dt 


d (d t ip*) 



5 e ip* 



d(p ih p) 



5* (PuP) + 



d A 



d{V lh S) 



S e {V.S) 



d n xdt . 



(A.9) 



For a fixed value of the index k and multi-index using k times Eq. (|A.8| ), we have: 



d C k 
o /^ NL x 5 e (Vj P) = E (-l) P A I r k V 1 



d /L 



*('0(2Xp). 



(A.10) 



where the coefficient v4 7 fe = n"=i ^'/(U m r !), S[ p =o] = Sp=o S Ip and the multi-index 
(ii, ■ • • , i n ), I p = (li, ■ ■ ■ , l n ) and I q = (mi, • • • , m n ) are related by i r = l r + m r . 
Using Eq. ( |A.10| ), Eq. ( |A.9| ) transforms to: 



5 e A 



dC L s ^ + d_ 
dip dt 



dC T 



d 

dxi 



d(d l ip) 



S e ip 



d(d t iP) 
dC 



Selp 



d_ 

dt 



dC T 



5 ip* + ® 
dip* e dt 



d(d t iP) 
dC 



S e ip + 



d 

dx,. 



dC L 

d(d t ip) 



d (d t ip* 



S e iP* 



d 


\ 1 


~ It 


d (d t ip*) 



6e1p 
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+ E E(-i) p <xR 

[fe=0] [p=0] I 



9 A 



5 A 



0(2>,S) 



5 t S\ \d n xdt . (A. 11) 



After inserting in Eq. ( |A.11| ) the expression of d C NL /dijj and d C NL /d ip* obtained from Eq. 
( |A.6| ) and its conjugate, we finally obtain: 



fc-i 



d(d t ^) 



+ E E(-i) p <^>, P 

[fe=l] [p=0] I 



9 (d t ip*) 



+ 



5 



dxi 



<9£„ TT 



p so J 



9 

<5 £ sUd n xcft . (A. 12) 



In presence of a symmetry the variation of the action must vanish and thus, from Eq. ( JA.12 ) 
after rearranging the terms, it follows the continuity equation: 

dQ 



dt 



V-JF = 



(A.13) 



with charge Q: 



9 (QeV" 



(A. 14) 



and the flux JF: 



T 

3 



d (djip* 



+ E E(-i) p <2\R 

[fe=0] [p=0] I. 



d JCL 



5 f S 



(A15) 



where T>. /fe a = dj with a = p, S, the coefficients = (ij + 1) A* /(m,j + 1) and 
= n — J2r=ij fio,m r ■ Recall that from the continuity equation the current is defined modulo 
the curl of an arbitrary function. This fact was taken into account in the expression of the 
current (|A.15|) . 

In the following we discuss two important cases. In the first, we suppose the system 
is [/(l)-invariant. Using the transformation ip — > tp exp(ie) where e is the infinitesimal 
generator, we have: 
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(A.16) 



5 € ijj = ieijj , 5 e ip* = —ieip* 

5 t p = , 5 e S = he . 

From Eqs. ( |A.14| ) and ( |A.15| ) we obtain the expression of the conserved density Q = p and 
the related current Ti = jf. 

du[ P , sy 

f U T 

[k=0] 



diS _ (-i) fc ^ 
— p+ E -7rdr v i 



(A.17) 



It is trivial to note that Eq. ( |A.13| ) is the continuity equation for the field ip, where the 
current ( |A.17| ) assumes a nonstandard expression, due to the presence of the imaginary part 



of the nonlinearity in the evolution equation. By taking into account the definition (|2.4j ) of 
the functional derivative, Eq. (|A.17|) can be written also in the form: 

U[p, S]d n xdt , (A. 18) 



Ji 



5 

p + 



m ' 5 (diS 

modulo a curl of an arbitrary function. Note that Eq. ( |A.18j ) can be obtained starting 
directly from the Eq. ( |2.7| ) after making the hypothesis that U[p, S] depends on the field S 
only through its spatial derivatives as required from the U(l) symmetry. 

In the second case, we discuss the energy-momentum tensor related to the space-time 
translations. Posing x^ — > x^ + e M we have: 

5 e ip = e M d^ip , 8 e if}* = e M d^ip* , 
5 e p = e M d^p , 5 e S = e M <9 M S , 

with p, — 0, • • • , 3 and do = d t . From Eqs. ( |A.14|) and (|A.15|) we obtain 



(A.19) 



oo 



2m 

h 



\ViP\ 2 + U{p, S] 



h 2 

(ditp* d t tp - diif> dtijj* 



(A.20) 
(A.21) 



T, 



i0 



2 m 



E E(-ir<^ 9 R 

[k=0] \p=0] 



d U[p, S] 



dtp + V Iv 



dU[p, S] 
d (V. T S 



dtS 



(A.22) 



T 



2 m 



(diip* djip + djif)* diip) + 5ijC 



E E(-i) p <?\K 

[fc=0] [p=0] 



OUjp, S] 
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OUjp, S] 



d 3 S 



(A.23) 



In Eqs. ( |A.20| ) and (|A.23|) we have taken into account the contribution due to the volume 
element. Note that the potential U[p, S] does not modify the expression of the momentum 
density Tqj which assumes the same form as in the linear theory. In contrast, U[p, S] 
changes the expression of the energy density Too, and more strongly the expression of the 
flux densities T ju . 



APPENDIX B. 

Theorem: If U[p, S] and F[p, S] are two smooth junctionals depending on the fields 
p, S and their spatial derivatives and satisfy the relation: 





dS 



U[p, S] = V • F[p, S] 



the functional U[p, S] takes the form: 



(B.l) 



U[p, S] = U[p, S] + V- G[p, S] 



(B.2) 



where U[p, S] depends on S only through its derivatives: dU/d S = 
Proof: Deriving Eq. (|B.2j ) with respect to 5* we obtain: 

d 



dV_ = dU_ d _ v 

OS dS dS ' dS ^ 

[k=0] I 



'» 'd(V.S) 



E 

[fc=0] L 



d(V Ik p) iV '* 'd{V Ik S)\ dS 

which coincides with Eq. ( |B.1|) for F = dG/d S. 
Alternatively, expanding the r.h.s. of Eq. ( |B.1|) : 



V ■ 



~dS 



dU 



E 



[fc=0] 



8 F OF 



(B.3) 



(B.4) 



div lkP ) - J d(v lk s) 

and integrating on the field S, after taking into account that p, S, T^ Ik P and T> 1 S are 
independent quantities, we have: 



u= E 

[k=o] 



OF d F 

^v^J-^-r + d^sy ; 



d(V Ih p) 



dtp.S) 



dS 
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S d ^ v ^Wr§-n) I F * dS+ S ^ s) a7?T^ / ^^ + (^ + (70^5 

[*=0] "l^PJ- 7 [*=!] °\ L 'i h ^) J 

V- / FdS + C- VS , (B.5) 



with Cj integration constants not depending on 5. Eq. ( |B.5| ) coincides with Eq. ( |B.2| ) for 
G = JFdS and TJ = C ■ VS. 

□ 
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